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Abstract 

>v 

' We investigate the evolution of rigid bodies in a viscous incompressible fluid. The flow 

is governed by the 2D Navier-Stokes equations, set in a bounded domain with Dirichlet 
boundary conditions. The boundaries of the solids and the domain have Holder regularity 
C 1 '™, < a < 1. First, we show the existence and uniqueness of strong solutions up to 
collision. A key ingredient is a BMO bound on the velocity gradient, which substitutes 
to the standard H 2 estimate for smoother domains. Then, we study the asymptotic 
Qh ' behaviour of one C 1,Q body falling over a flat surface. We show that collision is possible 

in finite time if and only if a < 1/2. 

1 Introduction 

To understand the dynamics of solid bodies immersed in a fluid is of primary physical 
interest, with regards to a wide range of phenomena such as sedimentation, filtration, 
or coagulation. For two-dimensional flows, under the assumption that the N bodies are 
T _j_ ■ rigid and homogeneous, and that the fluid is incompressible and viscous, one considers 

I/-) ' classically the following model: 

(N 

O n xeF(t). (1.1) 

qq ' div u = 0, 

O 



i) The velocity u and pressure p satisfy Navier-Stokes equations in the fluid domain F(t): 



ii) The N solid bodies are described by the closures S i (t) of connected bounded domains 
S l (t), 1 < i < N. They have rigid velocity fields 



u\t,x)=v l {t)+uj l (t){x-x l {t)) x , x£S l (t), l<i<N, (1.2) 

where v % and w l are the translation and angular velocities, whereas x 1 is the position 
of the center of mass. 

iii) The moving fluid and solid domains occupy a fixed bounded domain fl of R 2 , with 
Dirichlet boundary condition: 



F{t) = n\ufL 1 s t (t), u = o, xedn. (1.3) 
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iv) The fluid and solid systems are coupled by the continuity of the velocity, 



u = u\ xedS' l (t), (1.4) 
and the continuity of the stresses: 

r ( b \ r ^ ' 

J l uj l (t)= (x-x i ) ± -{^-pn)da+ (x-xY-p l f- 

JdSHt) \ on J 

The positive constants p, fi are the density and viscosity of the fluid. The positive con- 
stants m', J\ ~p l are the total mass, moment of inertia and density of the i-th solid. The 
source term / models an additional forcing (like gravity). The vector n at the boundary 
dU of an open set U refers as usual to the outward unit normal vector. 

Although natural, these equations exhibit some unexpected features, in both two and 
three dimensions. Hence, consider the case of one rigid body falling in a cavity (N = 1), 
under the action of gravity. It can be shown that if the boundaries of the body and 
the cavity are smooth, then no collision can occur in finite time. In other words, this 
system predicts that the kinetic energy of the body is strongly dissipated by the viscosity, 
resulting in no collision between the body and the boundary. This fact has been known 
from physicists for many years [H O H] , and was recently proved in one [32| and two 
dimensions [TUIISU] . 

This no-collision result is of course paradoxical. At the level of medium-sized objects, 
it goes against Archimede's law, and is clearly denied by common experiments. At a 
microscopic level, it also lacks relevance, as rebounds between particles are often involved. 
Many physics papers have been devoted to this paradox, trying to identify the flaw of the 
previous modelling. We refer to the articles [H [7] among many. Among the possible 
explanations, one of the most popular is roughness. Indeed, the no-collision result relies 
on the fact that the boundary of the solid structure is regular enough (namely C 1 ' 1 ). 
Small irregularities could then explain the occurence of collisions, see (26] [28] . Moreover, 
the effect of surface roughness in the dynamics of particles has been recently emphasized 
in experiments [HI H3J [33J • 

The aim of this paper is to study mathematically the roughness-induced effect on the 
collision process. Therefore, we consider Holder boundaries. Namely, we assume that 

dn e C 1 '", dS l eC ha , Vi, < a < 1. (1.6) 

We will first consider the well-posedness of system (|l.l[l - (|1.5|) . for such boundaries. We 
will establish existence and uniqueness of some strong solutions, up to collision. Our 
result extends previous results obtained for C 1,1 boundaries. Once this well-posedness 
is obtained, we will turn to the question of collision in finite time. We will consider the 
special case of one C 1,a rigid body, falling vertically over a horizontal flat surface. Losely, 
we will show the following: 

1. For a > 1/2, no collision can occur, and the strong solution exists for all time. 

2. For a < 1/2, one can find solutions for which collision occurs. 

This sharp criteria illustrates that roughness might be the reason for collision in fluid 
structure interaction, and the reason for the apparent paradox of the classical modelling. 

Before stating precisely the results, let us mention former mathematical studies. Fluid- 
solid interaction has been the subject of many papers, mostly devoted to the existence 
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theory for problem (|1.1[) - (|1.5[) . A key ingredient in many existence results is a weak 
formulation of the equations. Introducing the global quantities 



N 



v(t,x) := u(t,x)l F(t) (x) + ^u\t,x) l S i (t) (x), (1.7) 



N N 



p(t,x) := p F (t,x) + ^^(t.x) := ?lf(t)(aO + J2 ^ 1 S'(t) 0»0> ( L8 ) 

i=l i=l 

the conservations of global momentum, global mass, and bodies masses yield respectively: 
for all T > 0, for all <p € V, for all ip E D([0, T) x O), 



/ (pv ■ d t (fi + pv <E> v : D(lp) - 2{iD(v) : D(ip) + pf ■ ipjdxds + / po^o • ¥>(0) = 0, 

T 



(1.9) 



+ • w>) + po^(o) = o, 

[p%iI> + p 1 u • w) + / p^(o) = o. 



'0 J£l ' JU 

The space of test functions V is 

V = {^eP([0,r)xO), V-^ = 0, p\t)D{tp) = 0, Vi, VI < i < ivj. 

The divergence, rigidity inside the fluid, and no-slip condition read respectively: 

V-w = 0, p 4 D(v) = 0, 1 < i < N, v\ an = 0. (1.10) 

We refer to B. Desjardins and M. Esteban [10] for the derivation of these equations. 
Similarly to p, p l and w, the initial data po, Pq and vq are built upon the initial positions 
of the bodies Sq and the initial fluid and solid velocities uq, Vq, uJq. We will assume that 
there is no-contact initially, which means 

^n^ = 0, SlcQ, Vl<i,j<N, i^j. (1.11) 
Broadly speaking, previous studies deal with two kinds of solutions: weak and strong. 

Definition 1 A weak solution on (0,T), T > 0, is a family 

(S*(t),.F(f),i;), l<i<N, F(t) ^n\uf =1 SUj) 

such that 

i) S l (t) is a connected bounded domain, for all < t < T , for all 1 < i < N . 

ii) The scalar functions p, p l defined in (|1.8p and the vector field v satisfy 

(p,p J ) e L°°(0,T x a), v e L°°(0,T; L 2 (n)) n L 2 (0,T; ^(fi)). 

and equations (|1.9|) . 10[) . 

By classical results of R. Di perna and P.-L. Lions 12J on the transport equations Q1.9b .c). 
any (p, p l , v) satisfying ii) has the following additional regularity: 

P^eCMjL 1 ^)), 
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and the initial data is satisfied in this stronger sense. Moreover, any p 1 satisfying ([1.9b ) 
is the characteristic function of a measurable set: 

p % {t,x) = lgi( t )(:r), for a.e. t,x 

see [3S1 theorem 2.1, p. 23]. However, it is not clear that S l (t) should be open and con- 
nected, so that this constraint i) is added to the definition of a weak solution. Then, using 
the rigidity condition in (jl.lOp . one can deduce that v(t, •) is a rigid vector field on each 
S l (t), and by (|1.9b ). that >S l (i) = Ti-tSg, for a family of affine isometries IZt Lipschitz in t. 

The existence of global in time (T = +00) weak solutions was proved by E. Feireisl 
|14j and San Martin and coauthors [27] extending earlier studies "up to collision between 
solids" [TTJ [3TJ [55J 03 HE] ■ It holds in dimensions 2 and 3, with initial data satisfying 
and 

«o6i 2 (H), div Wo = 0, / € L 2 ((0,T); iJ _1 (fi)). 

Following the construction by E. Feireisl, no smoothness of the boundaries of the domain 
and the solids is necessary for the existence of weak solutions. However, the uniqueness 
of such solutions is unknown in general, even considering dimension 2 and pre-collisional 
times. After contact, it is known that uniqueness does not hold, as some entropy condition 
is missing to describe properly the post-collisional dynamics. This suggests to consider 
stronger solutions, namely 

Definition 2 A strong solution on (0, T) is a weak solution with the following additional 
regularity: 

v e L°° (Q,T;H^(n)) DL 2 (0,T; W hp (n)) for all finite p, d t v G L 2 (0,T; L 2 (n)). 

Our first result is the following: 

Theorem 1 (Well-posedness up to collision) 

Let v e H^(ft), p D{v a ) = 0, Vi, / G i 2 ((0,T); PF 1,00 (Q)), VT > 0. Assume 
(fl~TT[) . and 

dfleC 1 ^, a^eC 1 '", Vi, 0<a<l. 

Then, there exists a maximal € (0, 00] with a unique strong solution on (0, T) for all 
T < . Moreover, this strong solution exists up to the first collision, which means one 
of the following alternatives holds true: 

i) = 00, 6(t) > Vt. 

ii) T* < 00, 5(t) > Vf < T„, lim 5(t) = 0, 

where S(t) := min{ d(5*(i), S j {t)), d(S l (t),dn), l<i,j<N, i^j}. 

Note that by condition ([1.11)1 . and the Lipschitz dependance of IZt described above, S 
is positive at least for small times. Our theorem is an extension of results of B. Des- 
jardins and M. Esteban [TU], and T. Takahashi [31J . who proved respectively existence 
and uniqueness of strong solutions in the case a — I. See also [T7] for well-posedness 
under further technical assumptions on the solids. A key argument in these papers is the 
classical L 2 1— > H 2 regularity property for the inverse of the Stokes operator, which holds 
in C ' domains. In particular, one can show that 

f [ \S7 2 v{t,-)\ 2 < +00, 0<T<T*. (1.12) 

JO JF(t) 

In the case of general C 1,Q domains, this H 2 regularity result is still true away from 
the boundaries, and (p, v) still satisfies (|1.1[) in the strong sense, that is for almost every 
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Figure 1: Typical situation 



x, t. However Theorem 1 requires a control up to the boundary . We will show that the 
following BMO bound: 



substitutes to (|1.12|) . allowing for our well-posedness result. 

In a second part, we study if bodies can collide in finite time, that is if T* is finite or 
not. We consider one C < a solid that moves vertically near a flat horizontal surface undcr 
the action of gravity. More precisely let us denote S(t) — lZ t So the position of the solid 
at time t. We make the following assumptions: 

1. The source term is / = —ge 2 , with g > 0, e 2 = (0, 1). 

2. The solid moves along the axis X\ = 0, that is TZt is a vertical translation. 

3. The only possible collision points are on x\ = 0. 

4. Near x\ — 0, dil is flat and horizontal 

5. Near X\ = 0, the lower and upper parts of dS(t) are given by 



6. The solid is heavier than the fluid, i.e., P| s(t) > P\ F(t) - 

Note that if the initial configuration (Q, So, vq) is symmetric with respect to the xi-axis, 
then the unique strong solution will be symmetric for all < t < T*, and the solid will 
move along the vertical axis. Hence, there are plenty of configurations satisfying 1-6. A 
typical one is shown in figure [T] Our main result is the following: 

Theorem 2 (Link between collision and boundary regularity) 

For any strong solution satisfying 1-6, T* < oo if and only if a < 1/2. 

In physical terms, the theorem emphasizes the role of roughness in the collision scenario. 
Our result extends the results of M. Hillairet [20] and T.I. Hesla [19] in the case a = 1, 
for which it was shown that no collision occurs. Theorem [2] relies on the study of the 
stress J ds ^(fid n u — pn). When the boundary is regular, this stress diverges strongly as 
the distance to the boundary goes to zero. This mechanism prevents collision. When the 
regularity is weakened, the stress is also weakened, and contact may occur. The proof 
of the theorem involves the construction of appropriate test functions. In that respect, 




x 2 - = \xi\ 



l+a 



x 2 - x+{t) 



\ Xl \ 2 , 0<t<T„. 
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assumptions 2-5 are mostly technical, allowing to handle the computations. As can be seen 
from our proof, most of our arguments are local, and use only the weak bounds given by 
the conservation of energy. Hence, we believe that, as far as "real" (not grazing) collisions 
are concerned, the result might persist for more general domains and weak solutions. 
However, the source term must remain sufficiently integrable, as shown by an interesting 
example of Starovoitov [30j . Losely, Starovoitov exhibits an example of a weak solution, 
colliding in finite time, when Q and the solid are two spheres. But the corresponding 
source term satisfies only 

/ eL 2 (0,T; F _1 (0)), VT>0. 

The L 2 norm of f(t, •) diverges as 6 — > 0, t — > T*. This allows to compensate the divergence 
of the stress and to allow collision, even with regular boundaries. As shown by the first 
part of our theorem, this phenomenon is ruled out for more realistic forcing (such as 
gravity). 

The rest of the article is organized in three sections. Section 2 gathers regularity 
properties for the Stokes operator in C 1,a domains. Section 3 is devoted to the proof of 
Theorem Q] Section 4 contains the proof of Theorem [5] 

2 Regularity properties in C 1,a domains 

Existence and uniqueness of strong solutions have only been considered when solids have 
C ' boundaries. More precisely, a key argument in the papers of B. Desjardins and M. 
Esteban or T. Takahashi is the regularity estimate 

l|Vti|| H i(0) + \\p\\ H i(o)/m <C + IIsIIh^O)) (2.1) 

satisfied by the weak solution (u,p) of the Stokes system 

' -Ait + Vp = div F, xeO, 

div u = g, x e O, (2.2) 
u\do = 0, 

when O is a bounded domain with C > boundary. Such inequality is no longer valid when 
dO has weaker regularity. Only the interior estimate 

\\Vu\\ m{K) + \\p\\ H i iK y K <C(K) (\\F\\ H i {0) + \\g\\m(o)) (2.3) 

is satisfied, where K is any relatively compact open subset of O. As regards well-posedness 
issues, this interior bound is not sufficient. We will need a control up to the boundary, 
given by the following: 

Proposition 3 Assume that O has a C 1:Q boundary, < a < 1. Assume also that 

Fe L 2 (0)nBMO(0), g e L 2 {0) nBMO(O). 
Then, the weak solution (u,p) of (|2.2p satisfies 
||(Vu,p)||b M o(o) <C (\\(F,g)\\ 

bmo(O) + II {F, 9) IIl 2 (o)) ■ (2-4) 
We remind that BMO(O) is the set of functions / £ L x {0) such that 
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where the supremum is taken over all the open balls B of O, that is all the intersections 
of O with open disks. Note that the application 

IIZIIbmo(o) := su Pt4t / \f( x ) ~lB\ dx 

B \tS\ J B 

defines only a semi-norm, as it is invariant by the addition of constants. An easy remark 
is that BMO(O) is also characterized by 

sup inf f |-|- J \f(x) - m\ dx ^ < +00, 

where the infimum is taken over all real constants, providing an equivalent semi-norm. 
Note that / G BMO(O) if and only if / G BMO(R 2 ), where / is the extension of / by zero. 
Hence, standard results for the whole space apply directly to our setting. For instance, 
/ G BMO(O) belongs to L p (0) for any finite p, and 

sup f^j \f{x) - 7b\ P dx\ < +°°, 

this expression defining again a semi-norm which is equivalent to the previous one. We 
also remind the Sobolev imbedding in dimension 2: 

H\0) ^BMO(O), ||/||bmo(0) < C||/||jp(0). ( 2 - 5 ) 

which is simply deduced from Poincare inequality. Finally, we remind the interpolation 
inequality: for all 9 G (0, 1), for all 1 < p, q < +00 with (1 — 9)q = p 

||/IU»(0) < C 11/11 if fo) II/IIbmo(o) ( 2 - 6 ) 
We refer to [13, 24J for exhaustive study of the space BMO. 

Proof of the proposition. In the case of the whole space O = R 2 , the estimate (|2.4j) 
follows from the continuity of the Rlesz transform on BMO. In the case of a C 1,a bounded 
domain, it is connected to Holder theory for elliptic systems. Such theory has been of 
course widely considered, from various perspectives: see [TJ 1161 115] for some examples. 
Although a BMO estimate like (|2.4p may be part of the folklore of this domain, we could 
not find a proper reference for it. For the sake of completeness, we give here the main 
steps of (one possible) proof. The last step of the proof relies on ideas of Giaquinta et 
Modica, used to establish Holder estimates for the Stokes system with Neumann boundary 
condition [15] . 

We start with a simple remark, to be used implicitly throughout the sequel: any 
/ G L 2 (U), U open set, can be written / = div F, where F G i/ 1 (J7) satisfies 

II*1Ibmo(eo < C||^||ffi(to < C 11/11 L 2 (u) (2.7) 

This will allow to keep the source term in divergence form as we apply transformations 
to the Stokes system. 

Let (u,p) be the weak solution of (12. 2|) . where p is normalized such that J a p = 0. 
Standard energy estimates yield 

\\u\\hho) + blU'(O) < C {\\F\\ L 2 (0) + \\g\\ L 2 {a) ) . (2.8) 

Step 1 : Localization. Let O l <s O 1 , i = 1 . . . N a covering of O by open sets. Let ip l 
a smooth function with compact support in O l , such that tp' 1 = 1 on O l . The functions 

u 1 :— u, p 1 :— p 



7 



satisfy 

+ Vp' =div (^F) - FV# - 2(Vu)*V# - u + pV^ 1 

:=divF l , xeff, , 

(2.9) 

div u = gtp l + Vi> 1 ■ u := g\ x E 0\ 
u\ d0 i =0. 

By flU}, the L 2 n BMO norms of F i and g i are controlled by the I? n BMO norms of 
F and 3. Thus, we can restrict ourselves to a subdomain, that is establish (|2.4p with O l 
instead of O. 

Step 2 : Local coordinates. If % does not intersect the boundary of 0, the estimate 
follows from the interior regularity (|2.3[) . If O l intersects the boundary, we can assume 
with no loss of generality that it is a local chart: there exists a C l a diffcomorphism 

X : & ^ D(0, R), X {O l H dO) = (-R, R) x {0}, X (O* nO) = £> + (0, i?), 

where F> + (0, i?) is the upper half disk of radius R centered at the origin. We define new 
fields v, q, F', g' by the relations 

u\x) := v( X (x)), p\x) := q( X (x)), F*(x) = F'( X (x)), g l (x) = f^y- 
They satisfy 

-div (AVv) + div (Bq) = div (BF'), x G F> + (0, i?), 
B :Vv — g', x E D+(0,R), 

v \dD+{0,R) = 0. 

where 

^deW (Vx)tVx ' S ^dcrk)^- 

Note that A is uniformly elliptic over D(Q,R), and that A, S have C°' Q coefficients. As 
usual, (div M) 2 := and M : N = M 2j N tj for any 2x2 matrices M, N. 

Step 3 : Frozen coefficients. We write the previous system as 

-div (A(O)Vv) + div (B(0)q) = div (F), x G F>+(0, R), 

B(0):Vv = g, x G F> + (0, i?), (2.10) 

w lan+(o,i?) = 0. 
where 

F := BF' - (A(0) - A(x))V« + (5(0) - B(s))?, 5 := g' + (B(0) - B(x)) : Vv. 
Let us assume for a while that F E L 2 n BMO, g E L 2 P\ BMO, and that the estimate 

II (Vv,g) ||bmo(d+(o,jj)) < C (\\(F,g)\\ BMO (D+(o,R)) + ll(-P 5 5)II_l 2 (_d+(o,_r))) ■ (2- 11 ) 

holds. A simple scaling argument shows that the constant C can be chosen independently 
of the radius R. We now state the following a priori estimate: there exists a universal 
constant C, and s(R) going to zero with R such that 

\\(F,g)\\BUO(D+(o,H)) + \\(F,9)\\l2(d+(o,r)) < s(R) ||(Vt>,£?)|| B MO(z»+(o,Ji)) , 212 . 

+ C (\\{F',g')\\ B MO(D+{0,R)) + \\(F',g')\\L2(D+(0,R)) + || (V«, q)\\L3(D+(0,R))) 
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For the sake of brevity, we focus on the BMO bound, as the L 2 bound is straightforward. 
More precisely, we just show how to bound ||(A(0) — A(x))\7 v\\bmo, because the other 
terms composing F and g can be handled along the same lines. As emphasized at the 
beginning of the section, we need to control 

t b ■= t4t / I ( A ( x ) - Vv ( x ) ~ c \ dx 

\ B \ JB 

for any ball B of D + (0, R) and some constant vector c (possibly depending on B). Let r 
be the diameter of B and xo a point in B. We choose c = (A(x ) — A(0)) (Vw) B . We get 

Ib < ^ U \A{x) - A(x )\ (W) fl dx + J \A{x) - A(0)\ \Vv(x) - (W)^ dx\ 

< C [r a - 2 J \Vv{x)\dx + R a ||V«|| B MO(u+(o,fl))) 

< C (r a - 2+2 ^ \\Vv\\ LP{D+{0 , R)) + R a ||Vu||bmo(d+(o,ji») 

(2.13) 

for any finite conjugate exponents p,q, i.e. p 1 + q = 1. We choose q close enough to 1 
so that a — 2 + 2/q > 0. Together with the interpolation inequality (|2.6[) . we deduce that 

I B < C" (i? Q - 2+29 l|VHliV(o^))ll Vu llBMO(D + (o^)) + fl a HVt;||BMO(i 3+ (o,H))) 

< CiR 1 HVw||bmo(d+(o,7v,)) + C , 2||Vt;|| L 2 (£ ,+ (0!jj . ) ) 

for some universal positive constants 7, C%, C2. The estimate ()2.12|) follows. 

Note that estimates (|2.11j) and (|2.12|) yield the bound ()2.4|1 . Indeed, up to take smaller 
R, that is up to refine the covering of open sets O l , we can assume that s(R) < 1/(2C). 
Hence, combining (|2.12|) - (|2.1ip . we obtain 

h\ (Vw,g) ||BMO(23+(0,fl)) < \ ll(Vw,g)||L2 (£1 + (0 ^ )) 

+ CC'{\\(F',g')\\ 

BMO(D + (0,i?,)) + \\{F' 7 9'))\\l 2 (D+(0 : R))) ■ 

Then, it is well-known that L 2 , H 1 and BMO norms are preserved by C 1 diffeomorphisms. 
This allows to bound the right-hand side of the previous inequality: 

2 \\(^ v ,Q)\\l 2 (d+(o,r)) + CC (\\{F',g')\\ BMO{D+{0tR)) + ||(F',3'))IU 2 (c+(o,fl))) 

< C (\\ ul \\m(o*no) + \\p l \\L 2 (o*no) + \\( F \g l )\\L 2 (o*no) + II ( F \ S l )llBMO(0*ne>)) 

< C'(\\(F,g)\\ L 2 (0) + ||(F, fl )|| BMO (o)) 

where the last line involves the basic estimate (I2.8[) . As regards the left-hand side, we 
obtain the lower bound 

lb I ||BMO(CnO) < C IM|bMO(D+(0,.R)) 

and along the lines of (|2.13p 

l|Vu l || B MO(o*ne>) = ||VxVu(x(-))l|BMO(e?*ne>), 

< C (||Vw|| B MO(D+(0,i?)) + l|Vw|| L 2(D+(0,i?))) ■ 

This altogether implies ()2.4|) . 

We stress that (I2.11j) and (|2.12p are only a priori estimates: Vu,p, and therefore Vv, q 
are only supposed to be in L 2 , and not in BMO. Nevertheless, regularizing the coefficients 
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of A and B, establishing the same estimates for the regularized problem and passing to the 
limit allows to show that the weak solutions are indeed in BMO and that the inequality 
holds. As this regularization argument is very classical, we leave it to the reader. 

Step 4-' BMO estimate for the Stokes system in a half- disk. The final step of the proof 
is to derive the estimate (|2.1ip for the system (|2.10j) . By the reverse change of variables: 

x i ^ (Vx(O)*)^ 1 x 

we can assume that A(0) = B(0) = I2 is the identity matrix. By this linear mapping, 
the domain D + (0,R) turns into the intersection of a half-plane and an ellipse, say E + . 
As all the vector fields involved are compactly supported in E + , this Stokes system with 
Dirichlet boundary conditions still holds in any half-disk containing E + . As this system is 
rotationally invariant, we can furthermore assume the half-disk to be D + (0, R') for some 
large enough R'. Finally, as the estimate (|2.11[) is invariant by the dilations x 1— ► R'x, we 
can consider the case R' = 1. Eventually, we only have to establish the inequality 

||(Vu,p)||bmo(x>+(o,i)) < C (II(^3)IIbmo(d+(o,i)) + 9)\\l 2 (d+(o,i))) 
for the system 

-Au + X7p = div F, x£D + (0,l), 
div u — g, x £ D + (0, 1), 

u \dD+(0X) = 0. 



We remind that, if p is chosen such that 



p -- 

£>+(0,l) 



0, we already have the L 2 estimate 



II(Vu,p)||l2(b+(o,i)) ^ C \\{F,g)\\ L 2 {D+{0A)) . 

We shall rely on ideas of Giaquinta and Modica, who prove in article [TS] a Holder 
estimate for the Stokes equation with Neumann type boundary conditions. Let < p < 
R<2, and y in D+(0, 1). We will denote 



B(y,p) := %,p)nJ) + (0,l), f VtP := 
We decompose u 



1 



and (w,r) solves 



B{V,P)\ JB{y,p) 

v + 10, p = q + r, where (v 7 q) solves 

' -Av + \7q = div F, x £ B{y, R), 
div v = g - g y R , x € B(y, R), 

v \dB(y,R) = 0. 

-Aw + Vr = 0, xeB(y,R), 
div w = g y R , x £ B(y,R), 



f, V0<p<R. 



MdB(y,R) = u \dB(y,R)- 

We must first derive an estimate on v and q. We state without proof the well-known 
inequality (see [29] ) 

h-q y ,p\\LS(B(y, P )) < C\\Vv\\ L 2 {B{yiP)) , (2.14) 

where C does not depend on p by a simple scaling argument. Then, a standard energy 
estimate yields 



IB(y,R) 



F-F 



B{y,R) 



v,R) 



B(y,R) 



(9 - 9 v m) ■ (q - <ly, R ) 
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which combined with (|2.14p yields 

I|Vi;||i,3(b(j, ) ,r)) < C (\\F - F ytR \\ L 2 (B{y ^ R)) + \\g-g y>R \\&(B( v ,R))- (2-15) 
We now wish to obtain an estimate on w and r. As for q, the pressure r satisfies 

\\r -ry. p \\ L 2 {B{ytP)) < C\\Vw\\ L 2 (Bi y ip)) , (2.16) 
As regards w, we want to show the estimate 

p 2 

II Vw- &w)y lP U*{B{y,p)) < C-jp\\\7w~ {X7w) y R \\ L 2 {B{y , R)) (2.17) 

where C does not depend on p or R. At first, up to replace w by w — £2 ( ^^'v.r V 

which would still be zero at the flat part of the boundary dB(y, R) f] {x2 = 0}, and would 
still satisfy (|2.17p . we can assume that 

(d2Wi) yR = g VtR = 0. 

If p > R/2, inequality (|2.17| is trivially satisfied. If p < R/2, there are two cases. 

If B(y,R) C {a;2 > 0}, the ball B(y,R) does not intersect the boundary of D + (0, 1). 
We can use the interior estimate provided by Giaquinta and Modica in |15j : we can apply 
proposition 1.9, estimate (1.14) to the derivatives of w, which are still solutions of the 
Stokes equation, and this yields exactly (|2.17p . 

If B{y, R) n {x 2 = 0} ^ 0, we write 

\\Vw- (Vw) y p \\ L 2 {B(y ^ p)) < Cp || V 2 w\\ L 2 [B{ytP)) 

< C p 1 \\V 2 w\\ L ^ {B(y ^ p)) < Cp 2 \\V 2 w\\ L ^ (B{y ^ R/2)) 

< C{R)p 2 \\Vw\\ L 2 (B[ym 

Note that the first inequality is simply Poincare's inequality, whereas the last one stems 
from classical regularity results for the Stokes operator. Simple scaling considerations 
give the bound C{R) < C/R 2 for some constant C that does not depend on R. To prove 
(|2.17p . it is therefore enough to show that: for any solution w of the Stokes equation 

-Aw + Vp = 0, xeB(y,R), 
div w = 0, x(zB(y,R), 

satisfying moreover 

(d2W 1 ) yR = 0, w^Oon dB(y, R) n {x 2 = 0} (2.18) 

we have 

\\Vwh*my,R)) < C\\Vw- Ww) y , R h*(B(y,R)) if dB(y, R) n {x2 = 0} 7^ 0. (2.19) 
Again, the constant C in the r.h.s can be chosen independently of R. 

If inequality (|2.19p were not to be satisfied, one could find solutions w n satisfying 
(fZTTg]) . and such that 

\\Vw n \\ L 2 {B{ytR)) = 1, || Vw - {Vw n ) y R \\ L 2 {B(yM)) — -> 0. 
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From the first equality, up to a subsequence, w n —> w weakly in H 1 (B(y, R)). This implies 
the convergence of the averages (Vw n ) R — > (Vw) R . Moreover, by standard ellipticity 
properties of the Stokes operator, we have 

\\V 2 w n \\ L 2 {B) < C(R,B), VB^B(y,R) 

so that w n — ► w strongly in H 1 ^). Hence, we obtain, as in +00, 



V111 = (Vw) y R on B. 
From the second condition in (|2.18[l and the divergence-free condition, we get 



( d l W )y,R = °i ( 9 ^)yM = °' 



Moreover, by the first condition in (|2.i8[) , we also have {dzwijy R — 0. Hence, (Viy) R = 
0, and Vw = in any subset B relatively compact in B(y, R). Thus, Vio = on B(y, R) 
which contradicts the assumption that its L 2 norm is 1. 

This last argument leads to the desired inequality (|2.19p on (|2.17p . Combining (12.15|) 
and ()2.17p . we obtain 



||Vu - {Vu) y J\ L 2 {B{y>p)) < || Vw - (yw) v<p h*(B(v,p)) + ||Vu - {Vv) y J\ L 2( BM) 

/ p 2 _ 

C> \R2 H Vm _ ( Vu ) a ,_RlU 2 (B(?/,fl)) + W F ~ F y,R\\L*(B(y,R)) + ~ 9 y ,B,\\L*(.B(y,R)) 

< C" IIVU - (V^O^fllk^Bfe.R)) + IK^5)llBMO(D+(0,l)) i?: 

We use lemma 0.6 of [15] to conclude that 



< C 
< 



||Vu - (Vu)y iP \\L2(B(y : p)) < C||(F,p)|| BMO (£)+(0,l))P , 

which provides the BMO control of Vu. The BMO control of the pressure p then follows 
from 12~T3|) , 1(215)1 . This ends the proof. 

In the next section, we will use this proposition to show well-posedness of the PDE's 
system (ll.ll) - (|1.5|) . Before that, we state a regularity result of Sobolev type for the Stokes 
system in C 1:Q domains. It will allow to give a meaning in the trace sense to the stress 
tensor at the solid boundary (d n u — pn)\gs*- 

Proposition 4 Assume that O has a C l a boundary, < a < 1. Let s,t such that s < a 
and s < 2/t. Assume that 

F e L 2 (0) n W S ' T (0), g £ L 2 (0) n W S ' T (0). 

Then, the weak solution (u,p) of (|2.2[) satisfies 

\\(Vu,p)\\ W s,r (0 ) <C (\\(F,g)v\\ W s,r (0 ) + ||(F,.g)|| L2(0) ). (2.20) 

We remind that for all < s < 1, the fractional Sobolev space W S,T (0) is the set of 
measurable functions u such that 
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and this last expression makes it a Banach space. The assumption s < 2/t in the propo- 
sition ensures the continuous imbedding 

H\0)^W'< T {0), \\f\\ w ..r(o) < C||/|| K1(0) . (2.21) 
Similarly, the constraint s < a is such that C ' a (O) <-> W S ' T {0). 

Sketch of proof of the proposition. The proof of the Sobolev estimate (|2.20[) 
mimics the proof of the BMO estimate (|2.4[) , so that we only quote the few changes to be 
made. 

Steps 1 and 2 (localization and use of local coordinates) remain the same, up to the 
replacement of BMO by W s - T in every argument. 

In step 3, the only change is in the derivation of 

\\(F,9)\\w«--r(D+(o,R)) + \\(F,9)\\l 2 (d+(q,r)) < £(R)\\(Vv,q)\\w^(D+(o,R)) 
+ C (\\{F' ,g')\\ W s.T( D +^ R)) + \\(F',g r )\\ L 2( D +(p, R )) + ||(Vu, q)\\L*(D+(o,R))) 

which substitutes to (|2.12|) . Again, we just show how to bound ||(A(0) — A(x))Vv\\w s - T , 
as all other terms that compose F and g are treated in the same manner. We write 



\(A(0) - A(x) )Vv(x) - (A(0) - A(y))Vv(y)\ T 

D+(0,R)xD+(0,R) 



I D+(0,R)xD+(0,R) F V\ 

\V<y)\ Am -y V dxdy) 
D+(0,R)xD+(0,R) \x-y\ 2+ST 

< C I \\Vv\\ T Wa , T{D+{0tR)) sup \A{x) - A(Q)\ T + R T ^ \\Vv\\l T{D+(0>R)) ) 

< e'(R)\\Vvr w *, r(D+(0!R) y 

which allows to conclude as in the previous proof. 

Step 4, that is the W S,T estimate, < s < 1, for the Stokes equation in a half-disk, 
follows from a simple interpolation of similar inequalities for W°' T and W 1,T . 



3 Strong solutions 

This section is devoted to the proof of theorem [TJ Broadly, we shall prove existence 
and uniqueness of strong solutions as long as the distance between solid boundaries 5{t) 
satisfies 6(t) > So, where 5q > is arbitrary. The fact that strong solutions can not exist 
after collision will be discussed eventually. This altogether will of course imply the result. 
In what follows, constants will depend implicitly on 5q. 

We treat separately the existence and uniqueness parts. The existence result follows 
the lines of [TU] , whereas the uniqueness result is inspired by |H] . We thus rely substan- 
tially on these articles, and put the stress only on the changes due to our not so regular 
C 1,a boundaries. 

Our (refined) existence result reads: 

Proposition 5 (Existence of strong solutions) 

Let 5 > 0, v G i?o(0), PoD(v ) = 0, Vi, / G L 2 ((0,T) x Q), VT > 0. Assume 
(fTTT]) . and 

3(1gC m , 9Sj 6 C 1>a , Vi, 0<a<l. 
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Then there exists a strong solution on (0,T) for some T > 0. Moreover, one of the 
following alternatives holds true: 

i) One can take T arbitrarily large and S(t) > So for all t < T . 

ii) One can take T such that S{t) > So for all t < T and lim^T S(t) — So- 
In both cases, the strong solution has the additional regularity 



and 




where q is the corresponding pressure field, and 

F E (t) := {x e F(t) s.t. dist(x, dF{t)) > e}, e> 0. 

Proof of the proposition. Following [TO], we establish a priori estimates for a 
sufficiently smooth solution (v,q) on (0,T), s.t. S(t) > So for all t <T. 

We first take ip = v as a test function, which yields the standard energy inequality 

\\ v \\l°°(0,T;L 2 (Q)) + IMU a (0,T;ffJ(n)) < C(T) (\\vo\\ L 2( n) + ||/||l 2 ((0,T)x^)) (3-1) 

where C(T) is an increasing function of T. Then, we take (p = dtv as a test function, 
which yields 

/' / \d t v\ 2 +p, [ \D(v)(t)\ 2 <c( [ \D(vo)\ 2 + f [ \f\ 2 + f I \v-VvA (3.2) 
Jo Jn Jn \Jn Jo Jn Jo Jn J 

Note that the l.h.s. in ()3.ip . resp. ()3.2p controls the L°° n L 2 norm of v % , ui 1 , resp. the L 2 
norm of v l ,io l . We now use Stokes regularity to bound the last term in ()3.2p . 

The Navier-Stokes equation for the fluid part can be written 

[iAv — Vq = p(d t v + v ■ Vw — /) , div v = 0, x G F(t), 
v\ds*(t)=v i (t)+u i (t)(x-x i (t)) ± , v\ dnndF(t) = 0. 

As the solids and the cavity do not touch (8(t) > So), it is standard to build a solenoidal 
vector field w(t, •) € iP°(f2) such that 

w{t,-)\ SHt) = v l {t)+uo\t){x-x\t)) L , w(t,-)\ga = 0, 

with the estimate \\w(t, -)\\h s < C? a 5^-»(l' y *C*)l + I^WD- Then, the function u = v — w 
satisfies 

fj.Au — Vq = ~p~ (dtv + v ■ Vv — /) — fiAw, div u — 0, a; S -F(i), 

As F(t) is a C 1 '" open domain, we can apply the estimates of the previous section. If q 
is normalized so that j F u\ q(t, ■) = 0, we have, by (|2.7|) and propositions 13141 

||(Vu,?)(*)IU 2 (F(i)) + ll(VM,g)(0lli/i(^(t)) + II (Vu, <?)(*) || bmo(f(*)) 
+ ll(VM,g)(i)||iv>.-(F(i)) < C (l|9tw(i)|| L 2(F(t)) + II" • Vu(t)|| L 2 (F(t)) 

+ ||/(t)IU a (F(t)) + + |" < (t)|)). 
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We remind that this bound holds for all s, t such that s < a, and s < 2/r. Back to the 
original field v, and using the interpolation inequality (|2.6p . we get: for all finite r 

\\(Vv,q){t)\\ 2 Lr(m) + \\(Vv,q)(t)\\ 2 HHFB{t)) + \\(Vv,q)(t)\\l MO(m) 
+ \\(Vv,q)(t)\\ 2 ws , r{m) < C (\\d t v(t)\\ 2 L2{F{t)) + \\v ■ X7v(t)\\ 2 L2{F(t)) (3 4) 

+ Wmwhw)) + £(KWI 2 + KWI 2 ))- 

i 

By a time integration of (|3.4p from to t, and a linear combination with (13. ip and (13. 2p , 
we obtain 

/ (\\dtv(s)\\ 

+ \\(Vv,q){s)f ws , T{F{s) ^ds + ||V«(t)||| a( n) 

< C(T) (||«o|||x (n) + ||/||^ ((0> T)xn) + j[ll«-Vt;(*)||i a(n) d*) 

(3.5) 

where C (T) is an increasing function of T. To have a closed estimate, it remains to handle 
the nonlinear term. We split it into 



f / \Wv\ 2 = f [ \vVv(s)\ 2 ds + 

Jo J (I JO JF(s) i JO JS^s) 



The last term in the decomposition clearly satisfies 

rt r pt 

\vVv(s)\ 2 ds <CV |u;(s)| 4 + |u/(s)| 4 ds < C 




where C depends on ||t->o||l 2 and ||/||z, 2 ((o,T)xn) ■ The first term is bounded in the following 
way: 

/ \\v-Vv(s)\\ 2 L 2 (F r s)) ds < / \\v(s)\\ 2 LHFis)) \\Wv{s)\\ 2 LHF{s)) ds 
Jo Jo 

< C / ||u(s)|| L 2 (i , (s)) ||u(s)|| H i (F(s)) ||Vw(s)|| L 2 (F(s)) ||Vu(s)||BMO(F( ;i ))rfs 

Jo 

< C*||w|| L o=( ,t;L 2 (o)) / IKs)||Hi(n) l|Vu(s)|| B MO(F( s ))ds 

Jo 

< a + n I Hv«(«)||| MO(J r W )<i8 + c„ f \\Vv(s)\\i 2m 

Jo Jo 

(3-6) 

where 77 is arbitrary and C is an increasing function of ||vo||L 2 (fi) an( i H/llL 2 ((o,T)xfi) ■ 
Note that the second line is deduced from the use of Gagliardo-Nirenberg inequality and 
the interpolation inequality (|2.6p . Choosing 77 small enough, (|3.6p and ()3.5p imply that 

||V«(t)|U» (n) < C+ f ||V«(s)||i 3( n)||Vt;(a)||i a( n)ds. 
Jo 

using Gronwall lemma, and the fact that \\^v(s)\\ 2 ^ 1 , n -. is bounded through (|3.ip . we 
obtain 

l|Vv|| L =o( ,T,L 2 (n)) < C 
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where C is an increasing function of T, H^oH/f^fi) and ||/||l 2 ((o,t)xQ)- Using this bound 
in (13. 5p . we finally obtain: 




+ l|Vw(t)||i, (n) + \\(Vv,q)(t)\\ 2 H1(Fe(t)) + ||(Vf, 9 )(i)||| MO(F(t)) 



+ ll(Vu,g)(t)||w-».r(jr( t ))) dt + ||Vu||ioo( ,T,£ a (n)) < INIIiri(n), ||/IU 2 (o,Txn)) 

(3.7) 

These a priori estimates are as usual the key element in the construction of strong 
solutions, as it provides compactness for a sequence of approximate solutions. In the case 
of C 1,1 boundaries, the issue of building such approximate solutions and passing to the 
limit has been adressed in B. Desjardins and M. Esteban, as well as in many other studies. 
As it adapts straightforwardly to our case, we do not give further detail and refer to these 
papers. 

Let us stress that the W S,T regularity of (Vu, q) allows to define the stress tensor at 
the boundary (d n v — qn) \dF(t)- Indeed, taking indices s, r such that ts > 1 (together 
with the requirements s < a, t s < 2), one can define the traces of Vi> and q as elements 
of VK s-1 / T ' r (<9F(i)) for almost all t. Note also that the regularity properties 

v e L 2 (0, T; W M (n)), d t v e L 2 (0, T; L 2 (nj) 

of a strong solution v are enough to ensure that the right-hand side in (|3.3p belongs 
to L 2 ((0,T) x n). If 8(t) > S for all t < T, this automatically implies the L 2 (Hf oc ), 
L 2 (BMO) and L 2 (W S ' T ) bounds on (Vv, q) restricted to the fluid domain. This shows 
the last statement of the proposition, and concludes the existence part. 

We now turn to the uniqueness of strong solutions. Our result is 
Proposition 6 (Uniqueness of strong solutions) 

Let S > 0, v e H^(Q), Po D(v ) = 0, Vt, / G L 2 (0, T; W 1,00 (Q)), VT > 0. 
Assume p. lip , and 

dQ.eC 1 ' 01 , dSleC 1 ' 01 , Vi, 0<a<l. 
There is at most one strong solution on (0,T) such that S(t) > So for all t <T. 

Proof of the proposition. We follow closely the work of T. Takahashi related to 
C 1 ' 1 boundaries. We focus on changes due to our not so regular C 1,a domains. As in [3Tj . 
we just consider the case N = 1, / = 0, that is one solid S(t) immersed in the cavity fi, 
without forcing. To lighten notations, we also assume that the density p = 1 in the solid 
and the fluid domains. Minor changes allow to handle the general case. 

Step 1: Lagrangian coordinates. The first step in the analysis of uniqueness for this 
free surface problem is to get back to a fixed domain, by a change of variables of lagrangian 
type. Let vq e i? 1 (fi) and S(0) the initial velocity field and solid position. We will denote 
by h(t) the position of the center of mass of the solid at time t. We can always assume 
that h(0) = 0. Let (v, q) a strong solution on (0,T) such that 5(t) > So for all t < T. 

We consider the same change of variables as in (3U paragraph 4.1, pl504]: as S(t) > 5q, 
a solenoidal velocity field A(t,x) is defined such that 

A(t,x) = 0, for x in an 6q/4 neighborhood of dCl, 

A(t,x) = h(t) + uj(t)(x — h(t)) , for x in an <5 /4 neighborhood of S(t). 
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Then, one considers the flow 

X(t, •) : SI -+ n, -X(y, t) = A(t, X(t, y)), X(0, y) = y. 
which maps 5(0) to S(t) and -F(O) to F(t). More precisely, in a neighborhood of 5(0), 

X(t,y) = h(t) + R e[t) y, 9(t) = [ w(s)ds, Rg = (™» "j/) ■ 

Jo 

and near dVl, X(t,y) — y. Note that, as (v,q) is a strong solution, h,9 e H 2 (0,T). The 
mapping X inherits the regularity estimate 

\\dtX(t,-)\\ H s m < C.(\h®(t)\ + \6®(t)\), Vi = 0,l,2, V s € N. 

We then introduce the new functions 

u(t,y) := (Vyj^t.X^yJJw^.Jf^y)), p(t,») := q(t,X(t,y)). 

where Y := X^ 1 denotes the inverse of X with respect to the space variable, and as usual 

(vk)« = v; ; . 

Following [31, paragraph 4.2, pl507], equations (|l.l|) - (|1.5j) turn into 

' d t u + Mu + Nu - iiLu + Gp = 0, y e F(0), 
div u = 0, ye F(0), 

u( y ,t) = R-e(t)h(t) + "(t)y x , ye 5(0), 
mft(*) = / ([j,Vu- p)ndy, 

JdS(0) 

Ju(t)= ([iVu- p)n ■ y^ dy, 

JdS(0) 

plus the initial condition 

u\ t =o(y) = u (t,y) := v (t,X(t,y)). 

We refer to [31] for the exact expression of the various operators. In short, (dt + M)u 
corresponds to the original time derivative d t v, Nu corresponds to v-'Vv, Lu corresponds 
to Av, and Gp corresponds to Vp. An important point is that 

Nu = u- Vu, Lu = Am, Gp = Vp near dfl and 5(0). (3.8) 

Indeed, we have X(t,y) — y near <9f2, so that the change of variables is trivial near 
the boundary of the cavity. Similarly, VX(t,y) — Rgn) near 5(0). As Navier-Stokes 
equations are rotationally invariant, we get (|3 . 8[) . 

Step 2: Stokes like formulation. The operators above involve the flow X(t, •), which 
depends on the solution u itself: hence, they are nonlinear. But A(0, y) = y, which means 
that for small time, nonlincarities are expected to be small. We shall therefore treat these 
nonlinear perturbations as source terms. We introduce 

H{t) := / R_ eis) h(s)ds 
Jo 
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and write the system as: 

' d t u- fiAu + Vp = f - Mu-u-Vu, yeF(Q), 
div u = 0, y£ F(0), 

u(y,t) = ff(t) + w(t)y x , ye 5(6), 



mH(t) 
Jto(t) = 



(fiV« — p)ndy + w(t), 

dS(0) 

(fjiVu — p)n ■ y 1 - dy. 



(3.9) 



as(o) 



where 
and 



/ := - (TVu _ « . Vtt) + /i(£ - A)u - (G - V)p, 



«j(t) = muj(t)Ro it) h(ty 



Step 3: Uniqueness. The uniqueness of the strong solution will be established thanks 
to the formulation (|3.9|1 . Let (u 1 ,^ ), (u 2 ,g 2 ) two strong solutions on (0, T), T > 0, 
corresponding to the same initial velocity field vo E iJ 1 (fi) and same initial configuration 
5(0), F(0) = Cl \ 5(0). We remind that for the sake of brevity, we consider the force-free 
case. We assume that the boundaries dfl and <95(0) are C 1,a , and that 

S 1 (t)>S 0) S 2 (t)>S , Vt€[0,T). 

We can associate to v l the change of variable X 1 , the new functions u l , p 1 and so on. We 
shall prove that u 1 = u 2 on (0, T) x f2. The differences 



u — u 2 , 



H := H 1 - H 2 



uj x -uj 2 



satisfy with obvious notations 

' d t u - fiAu + \7p = f 1 - f + M 2 u 2 - M 1 !! 1 + u 2 ■ Vu 2 - u 1 ■ Vn 1 , y G F(0), 
div u = 0, ye F(0), 



u(y,t) = H(t) + u{t)y^, y G 5(0), 



mH(t) 
Juj{t) = 



(/iVw — p)ndy + w 1 ^) — w 2 (t), 

dS(Q) 

(/iVw — p)n • y 1 " ciy. 



SS(0) 



(3-10) 

with initial condition u|t=o = 0. Now, we can perform the exact same estimates as those 
performed earlier to show existence of strong solutions. In particular, we get (see estimate 

(\\dMt)\\l Hn) + ||V«(t)||| r(n) + UVu^lHi^o)) + ||Vu(t)||^ MO(F( o)) 

+ lb(*)llffl(F s (0))) dt + l|Vli|U«>(o,T;I/>(fi)) 



< RHS := C(T) H/ 1 -/ 2 ! 



li 2 ((o,T)xn) 



[u 1 ■ Vu 1 - u 2 ■ Vu 2 \\ 2 



L 2 ((0,T)xO) 



+ IIAfV 



M2u2 \\h((0,T)xQ) + \\ wl ~ W 2 \\ L i(0,T)) 



where e is any constant lower than So/ A, and C(T) an increasing function of T. As usual, 
the pressure p is normalized so that J F ^ p = 0. It remains to estimate the right-hand 
side. 
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By the remark (|3.8[) . f 1 — f 2 has a support F e which is compact in F(0). The L 2 (i? 1 ) 
bound on (Wu,p), which was true up to the boundary for C < domains, holds in F s . 
Moreover, the L 2 ((0,T) x £1) estimate on d t u and the L°°(0,T; L 2 (Q)) estimate on Vu 
also hold. Hence, the same bounds as those derived in [3TJ Corollary 6.16, pl523] apply: 

II/ 1 -/ 2 IIl 2 (o,txo ) < cr 1 / 10 (||(v u ,p)|| L2(0 , T;ff i (Fe(0)) 

+ \\dM\L*((0,T)xn) + I|Vu||l~(0,T;L2(O))^ 

where C denotes here and in the sequel an increasing function of T and || f^o II j? 1 (n) • 
As the L 2 bound on dtu is still available, we deduce as in [3TJ corollary 6.16], that 

Ik 1 - w 2 \\ L 2 [H 2 ) < CT 1/2 ||9 t u|| L 2 ((0iT)xn) . 

We remind that Mu = ui(t) u 1 - + dtY ■ Vit, . Therefore, 

||MV - M 2 u 2 || L 2 ((0 T)xn) < ||M 1 u|| i 2 ((0 T)xO) + \\Mu 2 \\ L 2 {(0T)xn) 

< C^||(w 1 ,/l 1 )|| i 2 (0 T) ||u|| L 2 (0:T;ff l ) + ||(CJ, /l)|| i 2 (0iT) ||u 2 || L 2 (0 T . ff l ) ^ 

< cr^llVuiixjojiHip)). 

Finally, we must control the quadratic term 

u 1 • Vu 1 - u 2 • Vu 2 = u ■ Vu 1 + u 2 ■ Vu. 
Like in previous computations, we get 

\\ u ■ Vul Hz,2((0,T)xF(0)) < C 'll u llL~(0,T;L4(O))ll Vul |U° c (i 2 (fi)) [ II Vu 1 (i) || BMO(f (0)) 

< c ^/r\\Vu\\ 2 LXI ^ T . H i^y 



using the L 2 (BMO) bound on Vu 1 . Similarly, 

||u 2 • Vu|||2(( 0iT)xF ( )) < C||u 2 || 2 00(i 4 (n)) ||VM|| L oo( i 2( 0) ) / ||Vu(t)|| B MO(F(o)) 



< C\fr [\\u\\ 2 Lx(H1{n)) + ||V?i|| L 2(BMO(F(0))) 

The L 2 bound in S(0) is straightforward, and we end up with 

\\u l ■ Vu 1 - u 2 ■ Vu 2 || L 2 (0 TxO) < CT 1/4 (\\u\\ 2 Loo(H i {n)) + ||Vu|| L 2 (BMO(F(0))) 

Eventually, these inequalities lead to 

RHS < CT 1 / 10 (||u|Uoo (H i ( n)) + ||fltu|U»(£ a (n)) 



+ l|Vu|| i 2 (BMO(F(0))) + ||(Vtt,p)||i2( ff l (FE)) j 



with C an increasing function of T and ||^o||ff 1 (n)- By reporting this bound in (|3.1ip . we 
deduce that there exists a small To such that u 1 — u 2 on [0,To]. Moreover, To depends 
only on ||t>o||jfi(n) (decreasing as ||fo||j? 1 (f2) increases). 
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Global uniqueness follows. Indeed, we know that u 1 and u 2 are in L°° ([0,T], H (d)) 
for all times T such that Vt < T, <5 X (i) , (5 2 (i) > 5a- Hence, up to consider a smaller To, 
we can apply the above local uniqueness result on [To, 2To], then [2To, 3To] and so on up 
to reach time T. This concludes the uniqueness part. 

So far, we have shown existence and uniqueness of strong solutions at least up to the 
first collision. Theorem [1] asserts more, namely that no strong solution can exist beyond 
the first collision time. This result can be deduced from [301 Theorems 3.1 and 3.2b]. 
Indeed, V. Starovoitov has shown the following: suppose that two C l a solids S n (t) and 
S l2 {t), resp. a C x < a solid S\t) and the C 1 '" cavity Q, collide for the first time at t = T. 
Denote for t < T, 

h iu i 2 (t) := dist(5 41 (t),S h (t))) resp. hi(t) := dist (S l (t), Si) 

and assume that 

u e T°°(0, T; L 2 (Q)) n L x (0, T; W x *(fi)). (3.12) 
Then, hi,hi u i 2 are lipschitz continuous on [0, T] and, for example, 

1 p + 1 1 



dhi , . 



for almost all f < T. 



< Chi(tf\\u(t)\\ w i, P{n) , (3 



1 + a p p' 



In particular, if T is the first collision time, and the strong solution exists beyond 
it, the regularity assumption ()3.12p is satisfied for arbitrary p. Taking p large enough, 
one can assume that [3 > 1. We note also that, by hypothesis (jl.lip . h(0) ^ 0. Then, 
by integration of the previous differential inequality, we obtain h(T) ^ 0, which yields a 
contradiction. 

We emphasize that for strong solutions, (|3.12[) holds a priori only for T < T». Thus, 
the argument of Starovoitov does not allow to conclude on the occurence of collision. 
In the next section, we will exhibit configurations for which collision occurs. For such 
examples, we have: 

IMIw^cn) = °°- 



4 The collision/no-collision result 

This section is devoted to the proof of Theorem[2l We consider the simplified configuration 
described at the end of the introduction, assumptions 1-6. In this framework, the position 
of the solid is characterized by h(t) := dist ((0, X- (t)), dft) . Later on, it will be convenient 
to use a parametrization by h, i.e. the translated domains 

S h := 5(0) + (h-h(0))e 2 , h e R. 

Of course, S^n) = S(t). By assumption 5, the boundary of S(t) is C ' near its "upper 
tip" (0,x+(t)), so that contact is impossible at this point, cf. [2D]- By assumption 6, 
gravity pushes S(t) downwards. So, we can even assume that 

inf dist ((0, dO) > 0. (4.1) 

te(o,T„) 



Thus, collision can occur if and only if lim t ^r, h(t) = 0. We will show that it is 
equivalent to a < 1/2. The proof is based on the use of a quasistationary velocity field w 
and quasistationary pressure field q. By quasistationary, we mean that for all t < T*, 

w(t,x) = Wh(t)(x), q(t,x) = q h (t){x) 



20 



for some stationary fields Wh(-), Ph(') defined on fl and parametrized by ft. > 0. Moreover, 
they will satisfy 

w G C^dp, T*); H\VL)), Aw(t,-) G L p (F{t)), p small enough, «e (0,2;), (4.2) 
div iu = in f2, w|s(t) = e 2 , w|asi = 0. (4.3) 

as well as 

qeC 1 ([0,T^),L 2 (Q,)), Vg(t,-) G L p {F(t)), p small enough, tG (0,T*). (4.4) 
In particular, we can use w as a test function in the variational formulation to get 



o Jn 



(pv ■ dfW + pv (g) v : D(w) — 2pD(v) : D(w) — pge2 ■ w) 



p(t)v(t) ■ w(t) - / p(0)v • w(0) 
n Jn 



Note that by 



/ pge 2 w = pge 2 • e 2 + / P5V(x h-> x 2 ) • w 
Jn Jsit) J Fit) 



= Psg\S(0)\ + Pf9 x 2 e 2 -n = (p s - pf) g\S(0)\ 

JdF(t) 

where p s := /o|s(t), Pf := p\p(t)- We also write 

2^ / D(v)(t) :D(w)(t) = h(t) f - qn J ■ e 2 - / (/iAw-V 9 )u 



9F(t) V ^ n / ^f(t) 

:= h(t)n(h) - / (/zAw - Vg) • v. 

Thus, the variational formulation yields 

AT(ft(t)) + ( ps - PF ) g \S(0)\ t = R(t) (4.5) 
where N is the antiderivative of n that vanishes at ft(0), and the remainder is 

R(t) ■= [ [ (pv ■ dtw + pv ® v : D(w)) + [ p{0)v Q ■ w{0) - [ p(t)v(t) ■ w(t) 
Jo Jn Jn Jn 

+ 11 Aw(s,-)- Vg(s, •)•«(*.•)• 

Jo JF(s) 

Theorem [5] will be deduced from the following proposition: 

Proposition 7 One can find Wh • i— > K 2 , % : fi i— * K, sizcft iftat satisfy ()4.2|) . 
dlH]) , P~4"|) . and sucft iftai 

i) For ft > small enough 

- c < n(h) < C h~ s \ /3=-^-, c,C>0. (4.6) 

I + a 

ii) For aZZ f < T», 

< C(H^olU-) (l + V*) - (4.7) 
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Before tackling the proof of this proposition, let us show how it implies Theorem [5J 
If a > 1/2, then j3 > 1. We get from (|4.5p and point ii) of the proposition: 
N(h(t)) > ( PF - ps)\S(0)\t -C(l + Vt) 
By point i), we also get for h small enough 

N(h) < -C\ln(h)\ 

In fact, one can take h 1- " instead of | ln(/i)| when (3 > 1, i.e. a > 1/2. Combining those 
inequalities, we deduce 

C\lnh(t)\ < (p s - PF )\S(0)\t +C(l + Vt) < +oo, Vi < T, 

which means that h does not go to zero in finite time. Hence, T* = +oo and there is no 
collision 

If a < 1/2, then (3 < 1, and n 6 L 1 . Thus, N is continuous. As h(t) is bounded, we 
deduce from (f¥T5|) : Vt < T, < +oo, 

-oo < inf N(h(t)) < (p F - p s )\S(0)\t + C(l + Vt). 

tG(0,T.) 

If T» = +oo, and > pj?, one can let i — > +oo in the previous inequality. As the r.h.s. 
goes to — oo in this limit, it yields a contradiction. Thus, < +oo. This ends the proof. 

The rest of the paper will be devoted to the proof of Proposition [71 
4.1 Construction of the test function 

We mimic the construction presented in article |20j for C ' boundaries. We want a 
function u>h{x) such that 

div w h = in fi, w h \ Sh = e 2 , w h \ dn = 0. (4.8) 

We always consider < h < hu := sup 0<t<T h(t), as no other value of h is involved in 
our problem. 

By a change of coordinates, we can assume (0,0) € dil, i.e. X-(t) = h(t). By 
assumption 5, there exists S > 0, such that 

VxedS h nD({0,h),2S), x 2 =~fh{x l ) := h+\ Xl \ 1+a , 

where as usual D(x, r) is the disk of center x and radius r. Moreover, by assumption 3 

and mm , 

S min := inf dist(dS h nD((0,h),5) c , dSl) > 0. 

0,h<h M 

To describe u>h away from the origin, we introduce a smooth function ip = ip(x), x £ R 2 
such that 

tp = 1 in a (5 m i„/2-neighborhood of S^^, tp = outside a <5 m i„-neighborhood of ^(o), 
We introduce another smooth function \ = x( x )> a; G K 2 , such that 
X = 1 in (-5, <5) 2 , X = outside {-25, 25f . 
Finally, we set Wh — X7 ± (xi<ph), with 
tp h = 1 in Sh, 

t 2 f 2ro \ 

<p h = {l-x(x))<p{x 1 ,x 2 -h + h(0)) + X (x) — ^ ( 3 - — jAr inQ\S h . 
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Figure 2: Geometry of the possible contact zone 



See figure |4~T1 to clarify the main notations. Note that ifh and therefore Wh are regular 
up to h — outside 

Q M := ^ n {|xi| < 6} n {x 2 < 7fc0i)}- 
Singularities at ft = correspond to the second term in the definition of (fh ■ 



It is straightforward that Wh satisfies (|4.8jl . As <ph involves the boundary function 
7^, the streamfunction X\iph has regularity C 2 a in the fluid domain. Moreover, wh is 
continuous across the solid boundary, so that it belongs to C°°((0, ft at); VT 1,0o (0)). In 
the fluid domain, its most singular second order derivatives behave like . We deduce 
that w(t,x) = Wht t \{x) satisfies (|4. 2[) . We postpone to the appendix the proof of the 
following estimates: 

Proposition 8 There exists < c < C such that 

\\wh\\L 2 {n) < C, 
c < ft^feT||V^|| L 2 (n) < C, (4-9) 
\\Vw h \\ L <x,(n\n hJi) < C, 



sup |7h(xi)| 3/2 / \Vw h (x 1 ,x 2 )\ 2 dx 
xie(-s,s) wo 




(4.10) 



P7fc0l) 

7/l (a; 1 ) 2 |0 A io h (a;)| a da: < C. 

1-5 JO 

Besides these estimates on Wh, the control of n(h) and R(t) shall involve quantities of the 
type 

/ {fJ.Aw h - \7q h )w 

JF(t) 

where w £ Hq(SY) is divergence free and satisfies w\gs h — e 2 - We prove in the appendix 
the following estimate 

Proposition 9 There exists a pressure field /iHg^E C°°(0, /im! C 1 (^)) such that for 
all divergence free w G H^fl) satisfying w\gs h = e 2 . 



/ (nAw h - \7q h )w 

JF(t) 



< C\\w\\ Hi{n) . 
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4.2 Proof of Proposition [7] 



Thanks to the estimates of the previous section, we are able to control the functions n(h) 
and the remainder term R(t). 

By integration by parts, 

n(h) := 2p I D(w h ) : D(w h ) + / (p,Aw h - Vq h )w h . 
Jn jf h 

By estimate (|4.9p b. we have 

ch-^rfKj < ||Vtflfc||£2(n) < Ch-^i, 
and Proposition [5] leads to 



(pAw h - Vq h )w h 



F(t) 



< C\\Vw h \\ L 2 (n) < C £ + ellVitfhHi: 



Combining these last two inequalities yields point i) of Proposition [7J 

To establish point ii), we need to control each term in the remainder. Still using bound 
we have 



/ {^Aw h(t) -\/q h[t) )v 



< C||w(*)llfl2(n)- 



Integration from to t and Cauchy-Schwarz inequality lead to 



f I (ma 

J JF(s) 

We also get 



w(s,-) -Vq(s,-))v(s,-) 



< C\\v\\ L 2(Q ft . HH Q,))Vl < C(\\u \\ L 2 (n) )Vt. 

(4.11) 



p(0)v o ■ w(0) - / p(t)v(t) ■ w(t) 
n Jn 



< C||u|| L oc (ot;i 2 (n)) SUP \\w h (t)\\ L 2 {n) 

he(0,h M ) 

(4.12) 



< C(||u ||z,3(n)) 



To deal with the term involving dtw, we shall use the following general bound: for any 
h e (0, h M ) and any (p,v) e L°°(Q) x Hq(Q) we have, for any w e Hq(Q) : 



pv ■ w 



< C\\p\\ L °°(n) ||Vv||i,2(n) \\w\\L*(si\n hiS ) 



<5 rjh(xi) \ I/ 2 

/ \ lh ( Xl )\ 2 \w(x)\ 2 dx) 
-s Jo J 

This is a simple consequence of Cauchy-Schwarz and Hardy inequalities, and its proof is 
therefore left to the reader. Note that dtw — h(t)dhW^ t y The previous formula yields 



/ I pv - d t w 


< C sup \h\ [ 


Jo Jn 


[0,T,) JO 



fS ph(s)(xi) \ 1 / 2 \ 

J J hh( s )(xi)\ 2 \d h w h(s) \ 2 dxj ) ds < C(\\u \\ L2m )Vt, (4.13) 



24 



where the last inequality involves (|4.10[) b. Finally, to deal with the nonlinear term, we use 
another general formula, namely: for any h € (0,/im) and any (p,v) <E L°°(Q) x Hq(CI) 
we have, for any w £ iJg(fi) : 



pv ®v : D{w) 



<C||p|U-(n) ||V«||i a(n) [\\D(w)\\ L ^ {n \n hs) 
+ sup It/^i)! 5 / |VtS(x)| 2 daj] 

xie(-6,5)\ wo 

This formula follows from Cauchy-Schwarz inequality together with a refined Poincare's 
inequality. We refer to lemma 12 in [20] for all necessary details. We infer from this bound 
and (|4~10|1 that 



pv®v : D{w) < C(|K|| L 2 (n) ) 



(4.14) 



o Jfl 



Gathering (|4.11|) to (|4.14[) gives the bound on R(t). 
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Appendix : Proofs of propositions [8] and [9] 

In this section, we estimate the rate of divergence of various Sobolev norms of Wh as h 
goes to 0. As explained in section 21 Wh is regular up to h = in Sh and fl \ (Sh U £lh,s)- 
Hence, there holds: 

l|Vt«h|U~(n\n ft ,,) < C, 
and the rate of divergence of Wh is the one of its restriction to Cl^g i.e. the one of 
V ± (xiiph) where 

Vh{x) = — -r-r 3 - 2 — - — - , Vie Q htS - 
lh{xi) V lh{Xl)J 

Proposition [5] is then a straightforward consequence of: 

Lemma 10 Given (a,p,q) € (0,oo) 3 , the quantity: 

\xi \ p dx\ 
_ 5 (h + \x 1 \i+")* 

behaves like 

(i) ch^" q , ifp+l<q(l + a), 

(ii) cln(/i), if p + 1 = q(l + a), 
(Hi) c, if p + 1 > q(l + a). 

when h goes to 0, with c a constant depending only on (a,p,q). 

The proof of this lemma as well as the induced bounds on u>h are direct adaptation of 
[201 Lemma 13]. 

It remains to build the pressure field q^ in order to prove proposition [HI For simplicity, 
we assume now /i = 1. With the same notations as in section [H we set: 

f Xl t 

qn{x) = d 2 i(xiiph(x)) + 12 / — - — r^dt, VxeO. 



7fc(zi) 



We stress that 



/•Xi 

q h (x) = d 2 i(x 1 ip h (x)) - / d 2 22{tiph(t,x 2 ))dt, Vxeflh.s- 

Jo 

As for Wh, this pressure field is smooth up to h = in the fluid domain outside &h,8- 
Consequently, the rate of divergence of Aw/, — Vqh is the one of its restriction to this 
latter domain. Standard computations lead to : 

Aw h (x) - Vq h {x) = „ , v , \x Vxeilhs- 
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We recall that V 2 Wh € L p (Qh,s) for p sufficiently small. As H 1 (Qh,8) C L r (Qh,5 ) f° r 
arbitrary r < oo, the integral to be estimated in proposition is well-defined. Up to a 
truncation (which leaves aside a term that is regular with respect to h), we can assume 
w = in (0, \ Sh) \ &h,6- A fortiori: 



/ (Aw h - Vqh) w= (Aw h - X7q h ) 



w. 



m\s h 

After an integration by parts, accounting for w\ds h = (0, w 2 ): 



(Aw h - Vq h ) ■ w 



dn(xi ip h )w 2 nida 



0S h 



9n(xi ip h )(2d 2 w 1 - diw 2 ) 



Thanks to Lemma HUl one can check that dn{x\ (fh) is bounded uniformly in h in L 2 (Q,h.s) 
then shows H^ii^/iHlsjq^ 5 )- Moreover, the boundary term reads 



dn{xi tp h ) w 2 nida 



dS h 



< \\w2\\L°°(dn h j) 



( 7 ^i)) 2 1 + K(*i)) s 



where |7^(a:i)| ^ c|a;i| Q . So, this boundary term is again uniformly bounded by Lemma 
ITU1 This ends the proof of proposition E3 
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